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Synthetic spin-orbit (SO) coupling, an important ingredient for quantum simulation of many 
exotic condensed matter physics, has recently attracted considerable attention. The static and 
dynamic properties of a SO coupled Bose-Einstein condensate (BEG) have been extensively studied 
in both theory and experiment. Here we numerically investigate the generation and propagation of 
a dynamical spin-density wave (SDW) in a SO coupled BEG using a fast moving Gaussian-shaped 
barrier. We find that the SDW wavelength is sensitive to the barrier’s velocity while varies slightly 
with the barrier’s peak potential or width. We qualitatively explain the generation of SDW by 
considering a rectangular barrier in a one dimensional system. Our results may motivate future 
experimental and theoretical investigations of rich dynamics in the SO coupled BEG induced by a 
moving barrier. 
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I. INTRODUCTION 

Spin-orbit (SO) coupling plays an important role for 
the emergence of many exotic quantum phenomena in 
condensed matter physics [1, 2]. In this context, the 
recent experimental realization of SO coupled neutral 
atoms provides an excellent platform for the quantum 
simulation of condensed matter phenomena because of 
the high controllability and free of disorder [3-6] of cold 
atoms. By dressing two atomic internal states through a 
pair of lasers, a Bose-Einstein condensate (BEG) with 
equal Rashba and Dresselhaus SO coupling has been 
achieved [3, 7-12]. The static and dynamic properties 
of such SO coupled BEG [13-23] have also been inves¬ 
tigated. Notable experimental progress in SO coupled 
BECs includes the observation of spin Hall effects [24] 
and Dicke-type phase transition [10], the study of col¬ 
lective excitations such as the dipole oscillation [8] and 
roton modes [25, 26], as well as the dynamical instabili¬ 
ties [27] in optical lattices, etc. Recently, the generation 
of another type of SO coupling, the spin and orbital- 
angular-momentum coupling, was also proposed [28-30]. 

Moving potential barriers have been used in the past 
for the study of the superfluidity of ultra-cold atomic 
gases. For instance, by stirring a small impenetrable bar¬ 
rier back and forth in a condensate, the evidence of the 
critical velocity for a superfluid was observed [31]. When 
a wider and penetrable barrier was swept through a con¬ 
densate at an intermediate velocity, the condensate is 
filled with dark solitons [32]. 

In this paper, we study the moving barrier induced dy¬ 
namics in a SO coupled BEG. We find that a fast mov¬ 
ing penetrable barrier may generate a dynamical spin- 
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density wave (SDW) in the wake of the barrier. Static 
SDW, which was proposed in solid state physics by Over- 
hauser [33, 34], has been widely studied in many different 
solid state materials such as chromium [35, 36]. Our gen¬ 
erated SDW in a SO coupled BEG is induced by the mov¬ 
ing barrier and vanishes when the SO coupling is turned 
off. The spatial periodic modulation of the spin density 
is not static, i.e., the local spin polarization oscillates in 
time periodically, and could last for a very long time. 

The paper is organized as follows. Section H describes 
the model of the SO coupled BEG. In Section HI, we 
study the dynamics induced by a suddenly turned-on sta¬ 
tionary barrier or a slowly moving barrier. Section IV 
includes the main results of the paper. We generate a 
dynamical SDW with a fast moving barrier, study its 
propagation, parameter dependence, and finally explain 
its mechanism using a simple one-dimensional (ID) sys¬ 
tem. Section V is the discussion. 


II. THEORETIGAL MODEL 

The SO coupled BEG is realized by shining two 
counter-propagating laser beams on cold atoms [3] . Two 
atomic internal states can be regarded as pseudo-spins 
I t) = 1^ = = 0) and | |) = |F = l^mp = —1) of 

^^Rb atoms. To stimulate the two-photon Raman tran¬ 
sitions, the two lasers are chosen to have a frequency 
difference comparable with the Zeeman splitting hwz 
between two spin states. The experimental configura¬ 
tion and level diagram are shown in Fig. 1. In our 
simulations, we consider a realistic elongated BEG with 
Y = 10^ atoms in a harmonic trap with trapping frequen¬ 
cies ujx,y,z = 27r X {20,120, 500} Hz. The strong confine¬ 
ment along the 2 ; and y directions reduces the dimension 
to quasi-ID. We use Er = as the energy unit, 

where kp the recoil momentum along the x direction 
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FIG. 1. (Color online), (a) Experimental setup for the ob¬ 
servation of barrier induced dynamics in a SO coupled BEG. 
The realistic barrier, created by another external laser beam, 
is characterized by the width W (along x direction), barrier 
peak potential 14, and barrier velocity v. (b) Energy level 
diagram of the SO coupled BEG. Two internal states are cou¬ 
pled by the two-photon Raman transition through a pair of 
counter-propagating Raman beams. 


(i.e., the SO coupling direction). 

A barrier, which can be created by the dipole potential 
of another laser beam [31, 32], is suddenly switched on in 
the BEG or is swept from the left to the right side with 
a velocity of v ranging from 1 to 80 fimjms. The barrier 
peak potential is 5 to 25 Er^ which is much larger than 
the chemical potential of the system. The width of the 
barrier is at the order of iim. The external potential bar¬ 
rier sweeping through the BEG is modelled as a Gaussian 
potential of the form: 


R(r,t) = 14e ^ ( 1 ) 

where Wx = W is the Gaussian barrier width along the 
SO coupling direction. much larger than the 

BEG widths in these two directions, xq and v are the 
initial position and velocity of the barrier potential, 14 is 
the peak potential of the barrier. 

The dynamics of the SO coupled BEG are governed by 
the Gross-Pitaevskii (GP) equation: 

~ [^*50 + Virap + V (r, t) + (2) 

where the single-particle Hamiltonian with SO coupling 
is given by 

Hso =KLCTzex) E-(Tz-\-—ax, (3) 

2m 2 2 

where {i = x^y^z) are the Pauli matrices, S is the 
detuning of the Raman transition and Q is the Raman 
coupling strength. The trapping potential is of the form 
Vtrap = muj^x^f2 -\-mujyy^/2 -i-muj^z^f2. For the sake of 
simplicity, we consider a 2D geometry in our calculations 
by integrating out the 2 :-dependent degree of freedom in 
the G-P equation 2, which is valid because the strong 
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FIG. 2. (Golor online). Excitation and propagation of the 
spin-dependent wave in a spin-orbit coupled BEG by sud¬ 
denly switching on a stationary barrier at the center of the 
trapped BEG at t = 0 ms. Golumn (a) is for spin up com¬ 
ponent, column (b) is for spin down component. 14 = lOEr, 
W = 0.5ym. Because of the presence of SO coupling, the 
propagation of sound demonstrates an anisotropic behavior, 
i.e., spin up only propagates to the left side while spin down 
only propagates to the right side. 


confinement along z direction restricts the BEG to the 
ground state of the harmonic trap along the z direction, 
yielding 

'>p3D{r,t) = i;2D{x,y,t) ^ (4) 


The interaction between atoms is determined by the 
mean-field Hamiltonian 


Hi 


0 


(5) 


where the reduced nonlinear coefficients for the 2D 
system are ’ 9n = = 

-^(co+c 2 ) ^ harmonic oscillator character- 

•^^' 4 ' maz 


istic length is az = y/h/muuz and the 5-wave scattering 
lengths are given by cq = 100.86ao, C 2 = —0.46ao (ao is 
the Bohr radius). 

In most cases the many-body interaction is weak, 
therefore many interesting physics can be well under¬ 
stood from the single particle band structure, which is 
either a double well type (for D < AEr) or a single well 
type (for D > 4Er) for (5 = 0. For ft < 0.2Er, BEG stays 
in both wells and the two dressed states interfere to form 
a stripe pattern; for 0.2F^^ < ft < BEG chooses 

either of the two wells as the true ground state, which 
is usually called plane wave phase or magnetized phase 


with a finite spin polarization \{(Jz)\ = 1 — (D/dE^)^ 
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FIG. 3. (Color online). A moving barrier with low velocity 
V — 1 jimlms pushes BEG to its right. Vh — 15Er, W — 0.5 
jim. Column (a) is for spin up component, column (b) is for 
spin down component. 


for Q > BEC condenses dit kx = 0 and the spin po¬ 
larization is zero \ {(Tz)\ =0. In our calculations, we focus 
on the latter case and take O = GE'r, S = OEr, where the 
generated SDW could be identified easily. 

One of the effects of the SO coupling is to change the 
sound of speed of the condensate. For a regula r BEC, 
the speed of sound is given by Vg = ^JUp/m^ where 
U = AttH^Ncq/ m is the nonlinear coefficient and p is 
the condensate density. When the BEC is dressed by the 
Raman lasers, the speed of sound is modified by changing 
the atomic mass to effective mass rrieff since the band 
structure of the system is modified, i.e., Vg = ^JUp/rrieff. 
For experimentally relevant parameters, the sound of 
speed is to the order of 1 pmjms. The sound of speed 
and the collective excitation spectrum has been measured 
in recent SO coupling experiments [25, 26]. In order to 
generate SDW, the velocity of the moving barrier should 
be much larger than the speed of sound. 


III. EFFECT OF A SUDDENLY TURN-ON 
STATIONARY BARRIER AND A SLOWLY 
MOVING BARRIER 

Before discussing the generation of SDW in the SO 
coupled BEC through a fast moving barrier, we consider 
two different limits: one is that a stationary barrier is 
suddenly switched on in the middle of the condensate, 
the other one is that a slowly moving barrier is swept 
through the BEC from the left side to the right side. 
Previously, suddenly switching on a barrier potential in 
the middle of the superfluid is usually used to measure 
the speed of sound for BECs or Fermi gases [37, 38]. 


In our calculations, the barrier potential is strong and 
therefore induces strong perturbations to the condensate 
which might be observed in experiments. Without SO 
coupling, the barrier excites two wave fronts propagating 
along both directions with the same speeds for the two 
spins. For SO coupled BEC, as shown in Fig. 2, the prop¬ 
agation of the two wave fronts propagate differently with 
respect to the direction of the spin. We see that only one 
clear wave front propagates to the left (right) side for spin 
up (down) with a speed around ^ 10 pm/ms for the cur¬ 
rent geometry and atom number. This anisotropic and 
spin-dependent propagation of the density perturbation 
is a direct consequence of the SO coupling. Note that 
except the propagation of the wave fronts, a series of 
density modulations are excited in the mean time, which 
are due to the suddenly switched on barrier induced per¬ 
turbations and occur for a regular BEC as well. 

When the barrier is slowly swept from the left side to 
the right side of the BEC, it is impenetrable for the con¬ 
densate because the barrier potential Vb = 15Er is much 
larger than the chemical potential ^ 1 of the con¬ 
densate. Therefore BEC is pushed in front of the barrier 
with the excitations of similar density modulations for 
the two spins (see Fig. 3). 

IV. SDW FROM A FAST MOVING BARRIER 

In this section, we focus on a fast moving barrier with 
a velocity larger than the speed of sound. Because of the 
increasing relative velocities between the barrier and the 
condensate, the barrier is now a penetrable potential for 
the atoms. For a certain parameter regime, the barrier 
induces a dynamical modulation of the densities of the 
two spins in the wake of the barrier, while it does not 
lead to any observable perturbations in its front. 

A. Generation and propagation of SDW 

Figure 4 shows the density distributions of two spin 
components when and after a moving barrier with a fast 
velocity is swept through the condensate. The velocity 
of the barrier = 50 pmjms is much larger than the 
speed of sound. We see that the density oscillations of the 
two spin components are out of phase, thus the barrier 
generates a SDW. The SDW could last a very long time 
and does not relax in the trap if all the parameters remain 
unchanged. 

To characterize the SDW, we calculate the wavelength 
A (the distance between two peaks for one spin compo- 
nent) and the contrast = \n'^rnax-K,in\l\Knax+Kxin\ 
near the center of the BEC and plot them as a function 
of the barrier’s peak potential 14, barrier’s velocity v and 
barrier’s width W. As shown in Fig. 5, the wavelength 
of the SDW is roughly a constant as a function of bar¬ 
rier height and width. However, the wavelength is almost 
proportional to the velocity of the barrier. 
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FIG. 4. (Color online). A fast moving barrier with a velocity 
of i’ = 50fim/ms penetrates the BEG and generates a dynam¬ 
ical SDW. Vb = 15Er, W = 0.5/im. Column (a) is for spin up 
component, column (b) is for spin down component. 



FIG. 5. (Color online). The wavelength A and the contrast C 
of the SDW generated in a fast moving barrier as a function 
of (a,b) the barrier peak potential Vb for v = 50fim/ms, W = 
0.5/i?n, (c,d) the barrier velocity v for Vb = 15Er, w = 0.5/i?n, 
(e,f) the width W ioT v — 50/im/ms, Vb — 


This is easy to understand because the SDW is not 
static. Each local spin polarization is fast oscillating as 
a function of time, as shown in Fig. 6 where we plot the 
local spin polarization {(7z){x = 0) in the middle of the 
BEG as a function of time. Before the barrier moves to 
X = 0, the density is a constant. Right after the barrier 
is swept through x = 0, the densities of the two spins 
start out-of-phase oscillations, therefore there is a spin 
polarization oscillation. The oscillation period is roughly 
a constant T for a certain barrier potential V^. Con¬ 



FIG. 6. (Color online). Local spin polarization {az){x = 0) 
as a function of time for Fig. 4. 


sidering that the barrier is swept through the BEG, the 
left and right side of the BEG is perturbed consecutively. 
When the barrier moves at a constant speed, the wave¬ 
length of the SDW should be proportional to the moving 
velocity if T does not depend on the velocity significantly, 
i.e., X = V X T. Note that this relation fails to apply for 
much larger velocities where the spin oscillation at dif¬ 
ferent local points may be generated at very short time 
and the above physical picture does not apply. In our nu¬ 
merical calculations, we verify that this oscillation curve 
remains the same when the interaction strength is var¬ 
ied in a large region, showing that the dynamical SDW 
is a phenomena governed by the single-particle physics. 
However, the phenomena changes when the interaction 
strength is strong enough such that the speed of sound 
of the condensate is comparable to the moving barrier 
velocity. As we have studied previously in Fig. 3, the 
barrier becomes impenetrable in this limit and pushes 
BEG to one side. 

In Fig. 7, we demonstrate the slow propagation of the 
SDW when the barrier is suddenly turned off after it 
moves to the center of the BEG. Since the dynamical 
SDW is actually a local spin polarization oscillation in¬ 
stead of travelling waves, the removal of the moving bar¬ 
rier at the center of the BEG should stop generating SDW 
on its right side. However, because of the superfluidity 
properties of the BEG and nonlinear interactions, the 
unperturbed neighboring atoms will be eventually per¬ 
turbed and we see that a new SDW with much smaller 
amplitude and velocity propagates to the right side of the 
BEG. The propagation of the new SDW is in the order 
of the speed of sound as expected. 


B. Mechanism of SDW generation 


In this section, we provide a qualitative theoretical 
understanding of the SDW generation by considering a 
rectangular potential sweeping through a ID SO coupled 
BEG system. The rectangular potential can be written 
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FIG. 7. (Color online). Propagation of SDW in a SO coupled 
BEG. When the barrier moves to the center of the BEG, we 
suddenly switch off the barrier and observe the slow propaga¬ 
tion of the SDW to the unperturbed region. All the parame¬ 
ters are the same as Eig. 4. 


as 


V{x, t) = H Mvt -Xo-W)- m{vt - Xo)], (6) 

where M{x) is Heaviside function: 

“<*)={? x> 2 ^ p) 

Similar as the Gaussian barrier, here Vb is the barrier 
potential and W is the barrier width. 

Figure 8 shows density distributions of the two spins 
at the moment when the rectangular barrier moves to a 
position around X = Sjnm for different barrier widths. 
We see that the width of the barrier changes the gener¬ 
ated SDW significantly. When L is integer times larger 
than the SDW wavelength L = nX, there are only n well 
developed complete SDW oscillations inside the barrier 
and the BEG that behind the barrier seems to be unper¬ 
turbed at all (Fig.8b,c,d). When L and A are incommen¬ 
surate, the SDW is generated in the wake of the barrier 
as shown in Fig. 8a. The wavelength of the SDW is also 
proportional to the velocity as we have demonstrated in 
Fig. 5(c) for the Gaussian-shaped barrier, while rarely de¬ 
pends on the barrier width and height. All these features 
of the dynamics can be well understood in the following 
way. 

The condensate is initially prepared at the ground 
state with D = 6Er while the band structure has only 
one single minimum at = 0, i.e., the quasi-momentum 
of the SO coupled BEG is zero. The initial wave function 
at some point x is 


The effect of a moving barrier could be explained using 
a simple single particle picture. In the laboratory frame, 
the real momentum of spin up and spin down components 
is = kx E kL = kL and k]^ = kx — kL = —kL respec¬ 
tively because = 0. An external barrier that moves 
along the SO coupling direction has different relative 
velocities for the two spins and therefore induces spin- 
dependent dynamics. Take the traveling external poten¬ 
tial as the frame of reference, the momentum are then 
~ —mu/fi + kL for spin up and k^B = —mv/h — kL 
for spin down. In the presence of the fast moving bar¬ 
rier, the velocities of the two spin components will be 
changed. From the conservation of the energy, we have 

the new velocities for the two spins: k'^^ = — k‘^ 

and k'^^ = — — k‘^ where ks = ^/2mVbJ^ = 2/cl 

since Now converting to the laboratory 

frame, the velocity is k'^ = k'^^ + mv/h for spin up, and 
k'^ = k'^Q-\-mv/hioY spin down. We find that for the large 
barrier velocity, ^ which means the quasi¬ 
momentum for this states is now k'^ = {k!^ k'^)/2 that 

agrees with the appearance of new momentum states in 
our GP simulations. Because of different group velocities 
of two spins in the presence of SO coupling, the moving 
barrier drives BEG to a new nonzero quasi-momentum 
states which is equivalent to an effective detuning Ag. 

Now the SDW related phenomena could be modelled 
as a quench dynamics where the effective detuning Ag 
is suddenly added to induce a coupling between the two 
new bands. Gonsider the effective Hamiltonian for a two- 
level system (ignoring other irrelevant constants): 

~ ( n/2 -Ae/2 ) ■ 


We denote ^ then the evolution of the 

local wave function for the point xq within the potential 
is given by: 


'0(^0, t) = 


2^2(1-^) 


2 - 


Me I inu 

-2 + A + 


• ( 10 ) 


Ignoring the normalized factor and the small terms of 
the order Ag, a straightforward calculation gives the spin 
polarization at this local point. 


{a^{xo,t)) = Ci^[cos{net) - 1], (11) 


where Ci > 0. 

It is quite clear that when the barrier moves to xq with 
a fast velocity, it perturbs the local condensate n^ixo^t) 
by coupling the two new bands and thus induces a local 
spin polarization oscillation {az{xo,t)). Due to the fact 
that the perturbation is applied from left to right, there 
is a relative phase between neighboring atoms. Therefore 
the spin polarization for a point on the left of xq is 


'0(x, t = 0) 



{az{x,t)) = C'i-^[cos(f2ei - ^^° ) - !]• (12) 


1 

-1 


(8) 
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FIG. 8. (Color online). SDW generation from a rectangular 
barrier in ID SO coupled BEG system for different barrier 
widths (a) W = 2jLim (b) W = 4.25/xm, (c) W = 8.5/im, 
(d) W = 12.75/im. Other parameter are 14 = 4Er, v = 
bO/j.m/ms. (e) The dynamics of the system is equivalent to 
the precession of a spin under the magnetic field. Without 
barrier, the system is in the eigenstate of a horizontal mag¬ 
netic field Bx. A moving barrier is equivalent to suddenly 
quenching on a z component magnetic field Bz. The spin 
then does precession about the axis determined by the total 
magnetic field. When the barrier passes the point, the spin 
could be in any possible orientations (with SDW appears) or 
returns to the eigenstate of Bx (no SDW appears). Other 
parameters are: Q = 6Er, 5 — OEr, Ux — 20Hz. 


We see from the above equation that the local spin po¬ 
larization within the rectangular potential is always neg¬ 
ative and the period of the oscillation is Tin = 27rfi/Oe. 
Therefore the wavelength is \in = 27rfin/Oe. 

Note that the point which has a separation of n\in with 
the right edge of the rectangular barrier has a vanishing 
spin polarization. If the left edge of the barrier coincides 
with these spin polarization vanishing points, i.e., W = 
nXin, then a striking effect occurs as we have seen from 
Fig. 8b-d: the condensate in the wake of the barrier seems 
to be unperturbed at all. That is because the potential 
has been removed for these points, where the dynamics 
is now governed by the original Hamiltonian with Ae = 
0. At the same time the wave function returns to its 
eigenstate with {az) = 0. 

If VF 7 ^ nXin^ then even though the governed Hamil¬ 


tonian returns to the original one, the wave function 
is not in its eigenstate and the coupling between the 
old two bands continues with the new Rabi frequency. 
Therefore in the wake of the barrier, we have the period 
Tout = > Tin and the wavelength Xout = "^Tout 

which is a little larger than Xin because Qe is slightly 
larger than Q. Similar calculation shows that the spin 
polarization could be positive or negative in one period, 
in agreement with the GP simulation (Fig. 8a). Further¬ 
more, when W = {2n + l)A^^/2, there is the largest spin 
polarization oscillation amplitude behind the barrier. 

The dynamics of the system are equivalent to the pre¬ 
cession of a spin under the magnetic field (Fig. 8e). With¬ 
out the barrier, the system is in the eigenstate of a hor¬ 
izontal magnetic field Bx and thus does not process. A 
moving barrier is equivalent to suddenly quenching on 
a 2 : component magnetic field Bz. The spin then pro¬ 
cesses about the axis determined by the total magnetic 
field. When the barrier passes through the point, the spin 
could be in any possible orientations (with SDW appears) 
or returns to the eigenstate of Bx (no SDW appears). 


V. DISCUSSION 

In our calculation, we have chosen Q = 6Er. Similar 
SDWs can also be generated for ft < AEr^ where the 
initial state is polarized and the population oscillation 
amplitudes of the two spins are now different. If the 
barrier is even faster than the speed used in this paper, 
then all atoms are perturbed almost at the same time, 
therefore the local spin polarization oscillation may be 
in phase. Since the delay is always present no matter 
how small it is, the pattern may look complicated. For a 
two-component BEG without SO coupling, a fast moving 
barrier does not induce any observable effects because the 
dynamics is governed by two uncoupled bands. 

The two-photon recoil momentum and recoil energy 
correspond to a velocity of Vr = hkLim = 4.\4[imjms 
and a kinetic energy of XE^. According to our simula¬ 
tions, to generate the SDW with a fast moving barrier 
(velocity n, peak potential I 4 ), we need to focus on the 
following parameter regime: 

KEb » Vb » Er, (13) 

where KE}y = ^mv‘^ is the corresponding kinetic energy 
of a particle with a relative speed of \ ± Vr — v\ ^ v 
(because v » Vr) respect to the barrier. 

In summary, we present a scheme to observe the 
generation and propagation of a SDW in a SO coupled 
BEG through a moving supersonic potential. The period 
of the SDW almost does not change with respect to the 
peak potential and width of the barrier. However, it is 
very sensitive to the velocity of the barrier. The essence 
of the SDW is due to the different group velocities of 
the two spin components in the presence of SO coupling. 
The SDW could last a long time in the trap without 
relaxation and therefore provides a good system to study 
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other complicated dynamics. For instance, by lowering 
the Raman coupling and changing the band structure, 
we may observe the opposite motion of the density 
modulations of the two spins and their relaxation in 
the presence of SO coupling. Furthermore, in other 
parameter regime (smaller barrier potential or velocity) 
or with a much narrower stirring barrier, it is possible 
to generate solitons or vortices. 
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